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With  allowance  for  galactic  perturbations  and  observational  error,  the  observed  distributions  of 
sizes  an  orientations  of  very  long  period  comets  are  consistent  with  a  uniform  distribution  of 
comets  within  the  Oort  Cloud.  ©  1985  Academic  Press,  Inc. 


Van  Flandern  (1977,  1978)  has  proposed 
that  the  very  long  period  “new”  comets 
originated  only  a  few  million  years  ago  in 
the  explosion  of  a  very  large  planet  in  the 
vicinity  of  what  is  now  the  asteroid  belt. 
The  primary  evidence  is  the  nonrandom 
distribution  of  both  sizes  and  orientations 
of  the  orbits  of  these  comets.  There  is  a 
sharp  peak  in  the  distribution  of  the  inverse 
semimajor  axes  (see  Fig.  la),  suggesting  to 
him  that  most  of  these  comets  passed 
through  their  previous  apparition  at  a 
common  epoch  (that  of  the  destruction 
of  the  planet).  The  nonrandom  distribution 
of  the  orientations  is  then  taken  to  re¬ 
flect  the  common  point  of  origin  of  these 
comets. 

Specifically,  there  is  a  suggestion  of  a 
nonuniformity  of  the  distribution  of  the  arcs 
traced  out  on  the  sky  by  these  comets,  pro¬ 
ducing  a  clustering  of  the  projected  inter¬ 
sections.  If  the  orbits  are  integrated  back¬ 
ward  one  revolution,  taking  into  account 
the  galactic  potential,  the  clustering  be¬ 
comes  more  pronounced  (see  the  figures  in 
Van  Flandern,  1977).  Van  Flandern  found 
the  point  of  maximum  concentration  (pre¬ 
sumed  to  be  the  point  of  origin)  to  be  at 
ecliptic  longitude  258°  and  “close  to  the 
ecliptic.”  Longitude  258°  and  latitude  0°  is. 
at  galactic  coordinates  /u  =  0°,  b11  =  +11°, 
i.e.,  close  to  the  galactic  center,  a  point 
noted  by  Van  Flandern  but  not  elaborated 
upon.  Ovenden  and  Byl  (1978)  have  con¬ 


firmed  this  clustering  at  the  previous  appa¬ 
rition,  although  they  have  not  interpreted  it 
to  indicate  the  origin  with  a  single  event. 

Recently,  Byl  (1983)  has  looked  more 
closely  at  the  effects  of  the  galactic  poten¬ 
tial,  and  he  has  pointed  out  that  one  result 
is  to  raise  the  perihelia  at  previous  and  sub¬ 
sequent  apparitions  to  the  point  that  pres¬ 
ently  observable  comets  would  be  out  of 
the  inner  Solar  System  at  those  times. 
Since  there  would  not  be  many  randomiz¬ 
ing  planetary  perturbations,  complete  equi- 
partition  of  the  cometary  orbital  energies 
would  not  be  possible.  However,  both  Byl 
and  Van  Flandern  seriously  underestimated 
the  degree  of  galactic  perturbations  and 
their  effects,  since  they  entirely  neglected 
that  part  of  the  potential  perpendicular  to 
the  galactic  plane. 

Consider  a  comet  moving  in  a  fixed  coor¬ 
dinate  system  with  the  x  axis  pointing  to¬ 
ward  the  galactic  center  at  the  present  ep¬ 
och,  the  y  axis  perpendicular  and  in  the 
plane,  and  the  z  axis  perpendicular  to  the 
plane.  Assume  that  the  comet  motion  is  suf¬ 
ficiently  limited  that  only  the  first  terms  in 
the  tidal  expansions  of  the  accelerations  are 
required.  Let  L  be  the  angle  from  the  x  axis 
to  the  radial  vector  toward  the  galactic  cen¬ 
ter  at  any  instant.  [L  is  180°  plus  the  galactic 
longitude  of  the  Sun  and  thus  has  a  period 
of  about  2  x  109  years.  Specifically,  L  — 
2.565  x  10-8  x  (year-1950.0)].  The  equa¬ 
tions  of  motion  can  then  be  written  as 
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x  —  CR[(1  +  cos(2L))x  +  sin(2L)y] 
y  =  CR[sin(2L)x  +  (1  -  cos(2L))y]  (1) 
z  =  Cz  ■  z. 

The  coefficients  have  the  following_approx- 
imate  values  at  the  solar  distance: 

CR  =  5  x  10“ 16  years-2 

Cz  =  -45  x  10-!6  years-2. 

[See  Mihalas  and  Routly  (1968)  for  a  com¬ 
plete  discussion  of  this  potential  and  the 
terms  in  it.]  The  third,  or  z,  equation  is  de¬ 
coupled  from  the  other  two  and  has  a  coeffi¬ 
cient  that,  while  an  order  of  magnitude 
greater  than  the  radial  one,  is  negative,  pro¬ 
ducing  a  high-frequency  oscillation  on  the 
radial  motion.  For  a  single  cometary  revo¬ 
lution,  this  is  the  dominant  term  and  must 
be  included. 

The  numerical  experiment  consisted  of 
taking  155  comet  orbits  with  initial  osculat¬ 
ing  parabolic  velocities  at  perihelion  and 
covering  the  entire  range  of  orbital  ele¬ 
ments  that  included  perihelia  within  the  in¬ 
ner  Solar  System,  and  numerically  integrat¬ 
ing  forward  one  revolution  those  that  were 
perturbed  into  long  period  eliptical  orbits, 
using  the  above  equations  of  motion.  All 
orbits  except  those  with  perihelion  direc¬ 
tions  close  to  the  galactic  radial  had  their 
perihelia  raised  to  well  outside  the  Solar 
System.  Those  comets  with  major  axes  and 
thus  virtually  all  of  their  motion  along  the 
line  from  the  galactic  center  experienced 
neither  tangential  nor  z  component  pertur¬ 
bations  of  any  significance,  thus  producing 
only  small  changes  in  their  orbits  (motion  in 
high-eccentricity  ellipses  is  essentially  rec¬ 
tilinear,  thus  making  the  initial  inclination 
irrelevant).  The  only  comets  that  therefore 
returned  to  where  they  could  be  observed 
again  were  those  that  had  closely  aligned 
lines  of  intersection  at  their  previous  appa¬ 
rition.  Independent  of  the  initial  distribu¬ 
tion,  presently  observable  comets  would 
show  a  clustering  of  intersections  toward 
the  galactic  center  at  their  previous  appari¬ 


tion,  making  this  a  rather  exotic  form  of 
observational  selection. 

As  a  secondary  effect,  the  lack  of  equi- 
partition  of  energy  means  that  the  distribu¬ 
tion  of  inverse  semimajor  axes  does  not 
have  to  be  uniform.  Indeed,  this  distribu¬ 
tion  must  almost  inevitably  be  peaked. 
Consider,  for  example,  the  case  in  which 
the  semimajor  axes  are  uniformly  distrib¬ 
uted  out  to  some  limiting  distance,  which 
will  be  assumed  to  be  a  sharp  cutoff  (which 
is  indeed  the  case,  due  to  both  the  galactic 
tides  and  stellar  encounters).  Define  x  to  be 
IQ6/ a,  where  a  is  the  semimajor  axis  in 
AUs.  If/(n)  is  the  distribution  in  semimajor 
axes  then,  in  this  example,  f(a)  ~  a2  out  to 
a  limiting  aa.  If  F{x)  is  the  corresponding 
distribution  in  x,  then  F(x)  ~f(a(x))(daldx), 
or  F(x)  ~  x-4  down  to  a  limiting  xo.  This 
distribution  increases  rapidly  toward  x0 
and,  after  convolution  with  a  normal  error 
distribution  for  the  observed  values  of  x, 
becomes  almost  indistinguishable  from  the 
result  by  assuming  all  true  x’s  have  the 
same  value.  A  test  was  run  by  assuming  all 


Fig.  1.  (a)  The  observed  distribution  of  IWa  for 
those  comets  clearly  classified  as  “new”;  (b)  the  simu¬ 
lated  distribution  in  the  same  parameter  from  the 
model  in  the  text.  The  values  of  l/aori8  were  taken  from 
Marsden  et  al.  (1978). 
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true  jc’s  come  from  a  single  parent  distribu¬ 
tion  that  can  be  characterized  by  a  single 
rms  error,  and  by  determining  by  trial-and- 
error  this  rms  error  and  the  x0  that  best  fit 
the  observed  distribution  of  x.  The  resulting 
predicted  distribution,  shown  in  Fig.  lb, 
results  from  an  rms  error  in  x  of  ±20  (106 
AU-1)  and  a  limiting  x0  of  25,  corresponding 
to  a  limiting  aphelion,  or  radius  of  the  Oort 
Cloud,  of  80,000  AU.  Both  of  these  esti¬ 
mates  are  quite  consistent  with  other,  quite 
independent  results. 

The  conclusion,  then,  is  that  the  ob¬ 
served  distribution  of  elements  of  very  long 
period  comets  can  tell  us  virtually  nothing 
about  the  complete  distribution  of  elements 
at  any  previous  epoch.  In  fact,  the  ob¬ 
served  distribution  is  not  inconsistent  with 
a  completely  uniform  distribution  of  comets 
within  an  Oort  Cloud  of  radius  80,000  AU. 
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